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ABSTRACT 

By means of numerical simulations, we investigate magnetized stellar winds of pre-main-sequence 
stars. In particular we analyze under which circumstances these stars will present elongated magnetic 
features (e.g., helmet streamers, slingshot prominences, etc). We focus on weak-lined T Tauri stars, 
as the presence of the tenuous accretion disk is not expected to have strong influence on the structure 
of the stellar wind. We show that the plasma-/? parameter (the ratio of thermal to magnetic energy 
densities) is a decisive factor in defining the magnetic configuration of the stellar wind. Using initial 
parameters within the observed range for these stars, we show that the coronal magnetic field config- 
uration can vary between a dipole-like configuration and a configuration with strong collimated polar 
lines and closed streamers at the equator (multi-component configuration for the magnetic field). We 
show that elongated magnetic features will only be present if the plasma- parameter at the coronal 
base is /^o <^ 1- Using our self-consistent 3D MHD model, we estimate for these stellar winds the 
time-scale of planet migration due to drag forces exerted by the stellar wind on a hot- Jupiter. In 
contrast to the findings of Lovelace et al. (2008), who estimated such time-scales using the Weber 
& Davis model, our model suggests that the stellar wind of these multi-component coronae are not 
expected to have significant influence on hot-Jupiters migration. Further simulations are necessary 
to investigate this result under more intense surface magnetic field strengths (~ 2-3 kG) and higher 
coronal base densities, as well as in a tilted stellar magnetosphere. 

Subject headings: MHD - magnetic fields - methods: numerical - planets and satellites: general - 
stars: pre-main sequence - stars: winds, outfiows 



1. INTRODUCTION 

T Tauri stars are recently formed low- mass stars (0.5 < 
M/Mq < 2), with spectral type ranging from F to M, 
and radius of no more than 3 — 4 Rq. They are the 
predecessors of solar-like stars that are still in the pre- 
main sequence (PMS) phase. Usually, T Tauri stars are 
classified in two categories: classical and weak-lined. A 
classical T Tauri star (CTTS) shows strong emission line 
activity, excess ultraviolet and infrared emissions, which 
are associated with accretion fiows and the presence of 
circumstellar disks. A weak-lined T Tauri star (WTTS), 
on the other hand, has a weaker emission line activity 
and shows little or no trace of accretion, implying that 
the accretion disk has been dissipated. For this reason, 
they are also known as naked T Tauri stars. 

WTTSs also show radio and X-ray fiare activity that 
is usually explained as a scaled-up magnetic activity 
similar to the one observed in the Sun. In the Sun, 
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highly detailed observations of the magnetic field, us- 
ing both in situ and remote measureme nts, have shown 
a ver y complex magnetic configuration (McComas e t al.l 
119951 : [Lin et al. 2005; Susino et al. 2008). Although in 
less detail, mag:n etic activity has also been observed in 
other stars (e.g., 'Petit et al."2005VPhan-Bao et al.'^ 
Reiners & Basri 2006; Wade et al. 2007; Johns- 
120071 : tValyavin et al. H2008l : iPonati et al. 2008). As ob- 
servational techniques are becoming more accurate, de- 
tection of a more complex magnetic field geometry in 
young stars has just started to b e accomplished, fo r 
example, the CTTSs y2129 Oph (|Donati et al.ll2007f ). 
BP Tau (IP onati et al.'^2QQg), and CV Cha & CR Cha 
(|Hussain e t al. 2009) where surface magnetic maps have 
been derived from spectropolarimetric data. Another 
system, V773 Tau A, a young binary system consisted 
of two WTTSs, is known to posses hi^h l evels of mag- 
netic activit y (IPhihips et all 119911 . Il996f ). Recently, 
Im assi et aP ()2008f ) observed two extended radio emis- 
sions that were interpreted as associated to the primary 
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and secondary stars of the system. They suggested 
that these features were produced by helmet streamer 
structures from each of the stars, that when overlapped 
at certain orbital phases, would result in reconnect ion 
producing the giant flares that were observed: radio 
emission would increase from a few mJy at apastron 
to more tha n 100 mJy at pe riastron. In another re- 
cent work, S kellv et alJ (|2008l ) reported the findings of 
slingshot prominences on TWA6, another WTT star. In 
their interpretation, such magnetic prominence extends 
for about 4 stellar radii, which is beyond the corotation 
radius located at 2.4 stellar radii. Magnetic field ex- 
trapolations from surface magnetograms have also shown 
that the geometry of t he lar^e-scale field of T Tauri stars 
may be very complex (iJardine et al.ll200"8l : iGregorv et al.l 
[2008). 

In the present work, we aim to investigate the mag- 
netic configuration of the corona of WTTSs using three- 
dimensional (3D) magnetohydrodynamics (MHD) nu- 
merical simulations. T Tauri stars are believed to 
be magnetized, and rotating several times faster than 
the Sun, hence it is important to consider a magneto- 
centrifugal wind. 

The first theoretical models to analyze the effects of 
both rotation and ma gnetic field on stellar wi nds are 
iWeber fc DavisI (|1967l . WD from now on) and iMestell 
(^968*), who investigated the coupling of a purely ra- 
dial stellar magnetic field to the wind using MHD ap- 
proach. Recently, theoretical and numerical efforts were 
made towards the study of more realistic magneto- 
centrifugal winds in diffe rent systems. For instance, 
IWashimi fc Shibatal (|1993l ) studied the influence of stel- 
lar rotation on the wind structure and acceleration of 
magnetized main-sequence solar-like stars investigating 
the effects of the centrifugal force. Holzwarth (2005), 
on the other hand, analyzed the impact of a purely ra- 
dial field on the wind. The stellar field, although ra- 
dial, was taken to have a non-uniform distribution in the 
stellar surface, as to account for the non-uniform sur- 
face magnetic fiux distribution observed in active stars. 
[Matt & Pudritzi (2008) aiming to explain the low rota- 
tion rates observed in CTTSs, studied the role of an ac- 
cretion powered stellar wind in the spinning down of a 
typical CTTS. 

All these works have confirmed that the structure of a 
magnetized stellar wind significantly d epends o n both the 
field intensity and topology. In Vidotto et al.l (2009), we 
performed an analysis of coronal winds in main-sequence 
solar-like stars. We showed that a decisive parameter in 
the acceleration of the magnetized wind of a solar-like 
star is the plasma-/? at the coronal base. Using 3D MHD 
self-consistent simulations, we explored how the change 
in the magnetic and thermal energy densities at the base 
of the wind modifies the magnetic configuration of the 
corona and the velocity profile of the wind. Rotation 
was not considered in that work. 

We now investigate the magnetic configuration of the 
winds of WTTSs. In our models we vary different phys- 
ical parameters of the systems, such as stellar rotation 
period, and magnetic field intensity and gas density at 
the base of the wind. In particular, we investigate the 
occurrence of helmet streamers or slingshot prominences, 
as well as their shape and extension. We perform self- 
consistent time-dependent 3D numerical simulations, in- 



cluding centrifugal, magnetic, gravitational, and thermal 
forces. The full set of MHD equations are solved to ob- 
tain the steady-state solution for the large-scale magnetic 
configuration of a TTauri star. 

The paper is organized as follows. In §2, we present the 
numerical scheme used, in §3 we describe our choice of 
parameters used in the simulations, and in §4, we present 
and discuss the results obtained. In §5, we analyze plan- 
etary migration due to drag forces exerted by the magne- 
tized stellar wind on "hot- Jupiters" , giant planets orbit- 
ing very close to the host star (< 0.1 AU). §6 is dedicated 
to conclusions and discussion. 

2. THE NUMERICAL MODEL 

To perform the simulations, we make use of the Block 
Adaptive Tree Solar-wind Roe Upwind Scheme (BATS- 
R-US), a 3D MHD numerical code developed at the 
Center for Sp ace Environment M odeling at University 
of Michigan (iPoweh et all IT999I ). BATS-R-US has a 
block-based computational domain, consisting of Carte- 
sian blocks of cells that can be adaptively refined for 
the region of interest. It has been used to simulate the 
heliosphere C ELoussev et al. 2003). the outer- helio sphere 
(jLinde et al.| [T998: Qph er et al.|[2QQl [2001. [20^ 
nal m ass ejections (Manchester et al."'2004'; 'Lusaz et al.j 
I2005D . the Earth's magnetosphere (Ridley et al. 2006) 
and the m agnetosphere of other planets (Toth e t al.l 
120041 : jHansen et al. 2005), among; others. In this work, 
we extend the model developed in lVidotto et al.l (|2009l ) 
to study the wind structure of WTTSs. 

BATS-R-US solves the ideal MHD equations, that in 
the conservative form are given by 
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where / is the identity matrix, p is the mass density, u the 
plasma velocity, B the magnetic field, p the gas pressure, 
g the gravitational acceleration due to the central body, 
and e is the total energy density given by 
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We consider ideal gas, so p = pksT / {prrip)^ where ks is 
the Boltzmann constant, T is the temperature, prrip is 
the mean mass of the particle, and 7 is the ratio of the 
specific heats. Due to the lack of knowledge of all the 
detailed processes that take place in a stellar wind, it is 
difficult to estimate all the mechanisms that modify the 
heat content of the wind (e. g., conduction, radiation, 
mechanical dissipation of energy that is transferred to 
the plasm a) . We adopt, therefore, an approach similar 
to Rousse v et al.l (^003) who considered that 7 is asso- 
ciated with "turbulent" internal degrees of freedom, in 
a way analogous to the Sun, where a significant amount 
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of energy is stored in the form of waves and turbulent 
fluctuations. 

For all the simulations performed, we adopted the same 
grid resolution. Initially, the simulation domain is re- 
fined in five levels. Other five refinement levels that are 
body-focused and focused on the equatorial plane (cur- 
rent sheet region) are applied next. Finally, an additional 
level is applied to the body. There are 9.1 x 10^ cells in 
the domain. The smallest cell size is 0.018 i?*, located 
around the central star, where is the stellar radius. 
The maximum cell size is 4.68 R^. The cell size near the 
current sheet is 0.036 R^. The grid is Cartesian and the 
star is placed at the origin. The axes x, y and z extend 
from —7 5 R^ to 75 R^ . More details of the grid can be 
found in'Vidotto et alj (|2009r ). 

The inner boundary of the system is the base of the 
wind at r = R^. Fixed boundary conditions were 
adopted at r = i?*. The outer boundary has outfiow 
conditions, i. e., a zero gradient is set to all the primary 
variables. 

The star is located at the center of the grid and has 
= 0.8 Mq and R^ = 2 Rq. The grid is initialized 
with a ID hydrodynamical (HD), thermal pressure driven 
wind for a fully ionized hydrogen plasma. This solution is 
dependent on the choice of the temperature at the base of 
the wind and on 7, and the only physical possible solution 
is the one that becomes supersonic when passing through 
the critical radius (Parker 1958). The initial temperature 
profile obeys the poly tropic relation T ex p^~^. Due to 
conservation of mass of a steady wind, we obtain the 
density profile from the radial velocity profile Ur{r). 

The star is considered to be rotating as a solid body 
with a period of rotation P. The rotation axis is along 
the z direction, parallel to the magnetic moment vec- 
tor. The simulations are initialized with a dipolar mag- 
netic field configuration described in spherical coordi- 
nates {r, ^, if} by 

B = (^cos^,isin^,0^ , (6) 

where Bq is the magnetic field intensity at the poles, r 
is the radial coordinate, is the colatitude, and (p is the 
azimuthal angle measured in the equatorial plane. As 
the magnetic energy density is latitude dependent and 
the thermal pressure is held constant at the base of the 
wind, the plasma- at the surface of the star is minimum 
at the pole. The MHD solution is evolved in time from 
the initial dipolar configuration for the magnetic field to 
a fully self-consistent non-dipolar solution, until steady 
state is achieved. 

3. THE CHOICE OF PARAMETERS 

In this section, we describe the choice of the parameters 
used in the simulations. They can be seen in Table [U 
where we present: the maximum surface magnetic field 
intensity Bq, the stellar period of rotation P, the ratio of 
specific heats 7, the density at the base of the corona po, 
and /3o, the surface value of the plasma- evaluated at the 
pole. We chose to simulate extreme values of the stellar 
rotation period, magnetic field intensity, and density, in 
order to probe the possible magnetic configurations for 
WTTSs. A brief explanation of our choices is presented 
next. 



TABLE 1 

The set of simulations. The columns represent, 
respectively: the name of the simulation, the 
surface magnetic field intensity at the pole, the 
period of rotation of the star, 7, the density at 
i?*, and the plasma-/3 evaluated at the pole. 



Name 


Bo (kC) 


P(d) 


7 


po (g cm 3) 


(3o 


A 


1 


3 


1.2 


1 X 10-11 


1/25 


B 


0.2 


3 


1.2 


1 X 10-11 


1 


C 


1 


3 


1.1 


1 X 10-11 


1/25 


D 


1 


10 


1.2 


1 X 10-11 


1/25 


E 


1 


3 


1.2 


2.4 X 10-10 


1 


F 


1 


3 


1.2 


2.4 X 10-12 


1/100 


G 


1 


0.5 


1.2 


2.4 X 10-12 


1/100 



Stellar Rotation: It is believed that, after the dissi- 
pation of the disk, a low-mass PMS star spins 
up as it contracts towards the zero age main se- 
quence. Such belief of a disk-regulated rotational 
evolution during the PMS phase is supported by 
observations of the bimodal distr i bution of rota- 
tion periods ( Bouvi er et al.l 1 19931 : iChoi fc HerbstI 
119961 : iRebulTet al.ll2QQ6h . In samples of PMS stars 
from the Orion Nebula Cluster and NGC2264, 
iCieza fc Baliberl ()2QQ7[ ) found that stars with disks 
possess a distribution of periods peaking around 
P ~ 8 d, while stars without disks presents the 
peak of the bimodal distribution at P ~ 2 d. To 
explore the extreme values of P, we chose periods 
ranging from 0.5 d to 10 d. 

Coronal Temperature: Over the last years, several obser- 
vations, such as from the lar^e X-ray surveys Chan- 
dra Orion Ultradeep Project (jGetman et al.l [20051 ) 
and XMM- Newton Ex tended Survey of the Tau- 
rus Molecular Cloud (G iidel et all 120071 ). showed 
that the coronal temperature of TTSs can exceed 
10^ K. It is believed that the bulk of the X-ray 
emission is coronal and it comes from the mate- 
rial confined insi de closed magnetic field structures 
(jPreibisch et al.l [2005). However, the outfiowing 
stellar wind is believed to be much cooler, as esti- 
mates presented by Joh ns-Krull fc Herczeg (20071) 
from UV spectra of the CT TS TW Hya and by 
IGomez de Castro fc VerdTI^ (|2007 ) for the CTTS 
RY Tau. In the present paper, our stellar wind 
models adopt a hot corona whose temperature at 
i?* is To = 1 X 10^ K. Future works exploring lat- 
itudinal dependence of coronal heating, such as a 
prescribed heating for loops, will be explored. 

Magnetic Field Intensity: High mean surface field 
strengths of up to a few kG have been observed in 
T Tauri stars (e.g [jQEns-Krulll[2007l : iDonati et all 
120071 : lYang et al.l I2008D . Although the magnetic 
field configuration at the surface seems complex, a 
dipolar component should dominate at larger dis- 
tances. A dipolar magnitude of few hundreds o f 
G has been detected ( Valenti fc Johns- Krulll[2QQl . 
In this work, we chose a magnetic field intensity in 
the range of 200 - 1000 G. 

Ratio of specific heats: When detailed heating and cool- 
ing mechanisms are unknown, the energy budget 
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of a stellar wind can be described by the use of 
a heating parameter 7 (Eq. [5]). Near the Sun, 
for instance, due to a large amount of heating 
(e.g., by turbulence, dissipation of MHD waves, 
etc.), th e solar plasma is considered as a gas with 
7-1 (jSteinolfso n fc HundhausenI 'l988). In T 
Tauri stars, the heating processes of the corona are 
poorly known. For simplicity, we parametrize the 
energy content of the wind of a WTTS by the use 
of 7 and we adopt 7 = 1.1 - 1.2. 

Density: The range of densities we adopted for the 
base of the stellar wind is 2.4 x 10~^^ - 2.4 x 
10~^^ g cm~^. This choice of density results 
in mass-loss rates ranging between — 10~^ and 
8 X 10-^ Mq yr-^ 

4. RESULTS: THE MAGNETIC FIELD CONFIGURATION 

In Fig. [TJ we present meridional plots for cases A, E, 
and F: the plots show contours of the poloidal velocity, 
the black lines are streamlines of magnetic fields, and 
the white lines represent the Alfven surface (where the 
poloidal velocity of the wind equals the local Alfven ve- 
locity). The first immediate conclusion is that the config- 
uration of the magnetic field lines is very different in each 
case. Case E (Fig. [1^) presents a dipolar-like configura- 
tion of the magnetic field lines. In case F (Fig. [it), the 
open-field lines are more collimated in the polar region 
and in the equatorial region, we observe the formation 
of streamers. In the intermediate case A (Fig. [TJd), we 
find a larger region of closed-field lines, coexisting with 
an open-field line region emanating from high latitudes 
of the star. What differs between these three cases is 
the plasma- i.e., the ratio of the thermal to magnetic 
energy densities. For cases E, A, F, /3o (evaluated at the 
pole) is respectiv ely 1, 1/25, 1/lQQ. 

As discussed in I Vidotto et al.l (|2QQ9l ) in the context of 
non-rotating main-sequence solar-like stars, the plasma- 

is a decisive parameter in the acceleration of the wind. 
When the magnetic energy density is of the same order 
as the thermal energy density, as in case E (Fig.[Tt), the 
magnetic configuration remains approximately that of a 
dipole. Otherwise, when the magnetic energy density 
dominates over the thermal energy density, the wind is 
accelerated reaching a situation where the flow ram pres- 
sure is able to distort the initial magnetic configuration, 
as in case F (Fig. [It). 

In addition to the acceleration of the flow, as a conse- 
quence of the Lorentz force, the wind velocity becomes 
more anisotropic with the decrease of Pq. Among the 
cases we analyzed, case E (and also B) presents veloc- 
ity profiles resembling spherical symmetry and in case F, 
where /3o = 1/100 (high magnetic energy density), the 
wind is visibly divided into a region of fast speeds at 
high latitudes and slow speeds at low latitudes. 

A comparison between cases A and F, where open-field 
lines are found, shows that: (i) case F shows higher colli- 
mation of lines towards the axis of rotation; (ii) the ratio 
of open to closed field line is larger in case F. As the 
magneto-centrifugal forces are responsible for the colli- 
mation of lines along the axis of rotation, it is expected 
that a wind more magnetically dominated would present 
higher levels of collimation (considering same stellar ro- 
tation rate). The reason why there are more open mag- 



netic field lines in case F is mainly due to the flow speed: 
the high wind velocities stretch and open the previously 
closed-field lines. As the polar wind is faster in case F, it 
is natural to expect a higher ratio of open to closed field 
line when compared to case A. 

These results also show that longer streamers will be 
present only if /^o <C 1. The presence of elongated mag- 
netic features has been recently considered to explain ob- 
servations of recurrent radio flares in binary s ystems com- 
posed of T Tauri stars (Massi et ajj i2008l : iSkellv e t al.l 
I2OO8I ). In these systems, the elongated magnetic fea- 
tures belonging to both stars could overlap in certain 
orbital phases, resulting in a modification of the mag- 
netic configuration (reconnection). This picture could 
account for giant radio flares, emitted when the binary 
system is near periastro n. Taking for example the system 
under consideration in Massi et al. (2008), the existence 
of elongated magnetic features indicates that /3o < 0.01. 
Adopting Bo - 2600 G and To - 1 MK, it results in 
noTo < 2 X 10^^ K cm-^ or po < 3.3 x 10"^^ g cm-^ 
where no is the number density at i?* and /i is the mean 
molecular weight of the wind particles. Therefore, the 
observations of radio flares can help us constrain the 
coronal density in these stars. 

PMS stars exhibit X-ray flares occuring in star-sized, 
small structures, often interpreted as sc aled- up solar-like 
coronal activities (' Preibisch et al.ll2005|). In addition, in - 
tense X-ray flares are also observed. iFavata "indl([2005h 
analyzed these intense X-ray flare decays to model physi- 
cal parameters of the flaring structure, flnding that these 
large flares take place in very long magn etic structures, 
extending out to several stellar radii. I Get man et al.l 
(2008) also argued that in fast rotating stars these pow- 
erful X-ray flares extend beyond the corotation radius, 
where centrifugal forces dominate over gravity. They 
conclude that the flares arise in traditional solar-type 
magnetic loops, whe re both footpoints are anch o red in 
the stellar surface. iJardine fc van BallegooijenI (|2005l ) 
developed a model for magnetic loops in rapidly rotat- 
ing main-sequence stars that explains the existence of 
magnetic prominences that can extend out to heights of 
2-5 i?*, above the corotation radius, in equilibrium with 
the open fleld region. In our simulations, we observed 
an extended, slender closed fleld region in case F. In this 
case, the extension of the streamer is assured up to an 
equatorial radius of ^ 13 R^. Above this point, the grid 
changes resolution, causing a local numerical reconnec- 
tion (in principle, the loop could extend farther out). At 
the equatorial plane, the corotation radius is given by 

m 

which for the TTS simulated here is r^o 
where Pd is the stellar rotation period measured in days. 
For P = 3 d, Tco — 4:R^. This shows that the elongated 
streamers, as obtained in case F, is stable to at least ~ 3 
times the corotation radius. 

Figure [2] shows the steady-state solutions for the re- 
maining of the simulations presented in Table [TJ We 
enumerate the following conclusions: {i) We show that 
only a change in the rotation of the star does not provide 
a signiflcant change in the poloidal velocity proflle of the 
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Fig. 1. — Meridional cuts of poloidal velociy profiles for: (a) case E (/3o = 1), (b) case A (^^o ■ 
lines are streamlines of magnetic field and the white lines represent the Alfven surface. 



1/25), (c) case F (^o = 1/100). The black 



system (compare cases A and D, F and G). It changes, 
however, how the magnetic field lines are twisted: the 
faster is the rotation of the star, the larger is the Lp- 
component of the magnetic field, (n) The heating pa- 
rameter 7 infiuences the thermal acceleration of the wind 
(compare cases A and C), being more accelerated in the 
case where the parametrized energy input is higher, i.e., 
for low 7. As the wind characteristics are different, the 
resultant interaction with the magnetic field is different: 
case A (high 7) presents a larger ratio of open to close 
lines than case C (case C has a larger closed-field line 
region), (m) Cases B and E present the same character- 
istics, as they possess the same (5^ (Vidotto et al.. .2QQ9). 
Furthermore, because /3o = 1, there is an equipartition of 
thermal and magnetic energy densities at the base of the 
corona, which results in lower velocities and field lines 
that are merely modified in relation to the initial con- 
figuration, resembling a dipole. {iv) Case G represents a 
star with a high rotation rate. Because of this, among the 
cases we analyzed, it is the case where the open field lines 
are more collimated towards the axis of rotation. Com- 
paring it to case F, where the rotation rate is 6 times 



smaller, it can also be seen that the higher the rotation 
rate is, the less elongated are the closed field lines. 

The Alfven surface plays an essencial role in the de- 
termination of the magnetic field configuration. In 
magneto-centrifugal winds, both B^jBj. and B^jBe are 
expected to be small inside the Alfven surface. Beyond 
the Alfven surface, the r and d components of the mag- 
netic field decay faster than the component, causing 
the field lines to be twisted. Therefore, helmet streamers 
anchored on rotating stars are not radial. In Fig. [3l we 
plot a 3D view of the region close to the star. As can 
be seen, the picture of approximately radial streamers, as 
observed in the Sun or in slowly rotating stars, is not any- 
more verified in WTTSs with high rotation rates: both 
open and closed lines become twisted due to the stellar 
r otation. 

Andre et al. (1988) presented an analytical expression 
for the calculation of the Alfven surface, considering that 
(i) the wind should have spherical mass-loss; (ii) by the 
time the wind is crossing the Alfven surface, it has al- 
ready achieved the terminal velocity; (iii) the field is 
dipolar all the way to the Alfven surface. As can be 





Fig. 2. — Same as Fig. ^for the remaining cases. 



seen from our results, when the self-consistent interac- 
tion of the field lines and the wind is taken into account, 
t hese hypothe s es ma y not be verified any longer. 

iMassi etld] ()2QQ8f ) computed the Alfven radius of the 
wind of each of the stars in the binary system V 773 
Tau A, whe r e the y observed helmet streamers, using 
I Andre et alJ (LlDSS^s analytical expression. They ruled 
out the twisting of magnetic lines of the helmet streamer, 
because in their estimate they found a large Alfven sur- 
face with radius of ^ 60 R^. According to our model, the 
Alfven surface of each individual stellar wind is located 
much closer to the star, and is probably even smaller 
than the minimum separation of the stars in the binary 
system V 773 Tau A (~ 30 i?*), suggesting that the in- 
teracting streamers are probably twisted, as illustrated 
in Fig. [3 

5. STELLAR WIND EFFECTS ON PLANET MIGRATION 

In the previous section, we investigated the magnetic 
configuration of the wind of WTTSs. This configuration 
arises naturally from the interplay between the outfiow 
and the field, i.e, the magnetic field is not restricted in 
our simulations. Several models assume a fixed magnetic 



field topology, as is the case of the Weber & Davis (WD) 
model. In this section, we investigate the effects our more 
realistic wind model will have on the migration of a giant 
planet. 

According to the current favored planetary forma- 
tion theory, the core of a giant planet forms far away 
from the central star, beyond the snow line (> 5 AU, 
where the disk temperature falls below the condensa- 
tion t emperature of water) by coalescence of planetesi- 
mals (iPollack et al.|[T996,) . When the solid core is suf- 
ficiently massive, it captures the surrounding gas of 
the acc retion disk to form the at mosphere of a giant 
planet (|Papaloizou fc Nelsonll2QQ5l ). However, observa- 
tions show the existence of planets orbiting solar-like 
stars at very close distances (< 0.1 AU). For this rea- 
son, these giant planets are often referred as hot-Jupiters. 
A possible explanation for this inconsistence is that the 
giant planet may have formed several AU from the cen- 
tral star, and later migrated inward. The interaction 
of the protoplanet with the disk wherein it was formed 
i s a probable p rocess that may drive this migration 
(|Lin et al.lll996f ). 

The inward migration is expected to cease when the 
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the planet's angular momentum leads to 



Fig. 3. — Three-dimensional view of the inner most region of case 
G. The hnes plotted were selected as to outline the twisting of both 
closed and open magnetic lines anchored on the surface of the star. 



planet is in the inner magnetospheric cavity, a cleared 
region between the stella r surface and the inner re- 
gion of the accretion disk (|Lin et alj|1996l : iTrilling et alJ 
[TqQSI: iPap aloizou 2007V In re cent numerical simulations, 
iRomanova fc Lovelacg (|2006f ) found that for a stellar 
dipolar field that is not highly misaligned with the spin 
axis of the star neither almost aligned so that 3D instabil- 
ities cannot fill the magnetospheric cavity with matter, 
the formation of the cleared region is assured. As in the 
cavity there is insufficient external material to remove 
the planet's angular momentum, the giant planet radial 
motions halt; the planet cannot grow further. In addi- 
tion to this mechanism, others have been proposed to 
explain why the p lanet migration is inhibited all the way 
to the star (iLin e t al. 1996; Fleck 2008). 

In this section, in particular, we will investigate the 
action of magnetic torques from the stellar wind acting 
on the planet. 

Considering that a T Tauri star has magnetic fields and 
rotation rates that are several times larger than in the 
Sun, Lovelace et al. (2008, LRB08 from now on) stud- 
ied the torque that a magneto-centrifugally driven wind 
would cause on a close-in giant planet. They adopted 
the model of WD for the stellar wind and evaluated the 
change of the planet's angular momentum due to an az- 
imuthal drag force exerted by the wind on the planet. 
Considering a solar- mass star with 2 Rq, P ~ 3 - 5 d, 
and Bq ^ 1 - 3 kC, they found a time-scale of 2 - 20 Myr 
for a giant planet like Jupiter, orbiting the equatorial 
plane of th e star, to have c onsiderable radial motions. 

FollowinglLovel ace et al.l ((2008), we estimate this time- 
scale by adopting a more realistic stellar wind, product 
of our simulations. Consider a planet of mass Mp, or- 
biting very close to the star at a distance Vp. The an- 
gular momentum of the planet is then Lp = MpVx^p^ 
where vk = {GM^/vpY^'^ is the azimuthal velocity of 
the planet, assumed Keplerian. Therefore, a change in 



dt 



lM,„.f^4M,„i, (8) 



where is the time-scale for the st ellar wind dra g 
change significantly the planet's orbit ( Lin et all 1 19961 ). 
Another way to express the torque (Eq. [8]) is to multiply 
the force the stellar wi nd will exe rt on the planet by the 



distance Vp ([Lovelace et al.|[2008f ) 



dt 



(-Pram^eff )^p 9 



(9) 



where Pram is the ram pressure the wind exerts on the 
planet 

1 . .0 1 



(10) 



and Aeff is the cross-section of the planet that intercepts 
the stellar wind. In order to take into account the ef- 
fects that the magnetic field of th e planet could have in 
the drag force of the stellar wind, iLovelace et al.l (|2008l ) 
defined an effective cross-section of the planet as 



eff 



max 



nRi 



1/3 



(11) 



where Bp is the surface magnetic field intensity at the 
pole of the planet and Rp is its radius. The planet is 
assumed to have a dipolar-field configuration, and plan- 
etary mass loss is not considered. From Eqs. ([8j) and ([9]), 
we find 

^^^ijv^. (12) 

2 Pram^eff 

A large Aeff means that the cross-section of the planet 
intercepting the wind is large, and thus, a more efficient 
drag mechanism (i.e., low r^) is expected. If the wind 
ram pressure is large, the dragging mechanism is also 
expected to be important. However, a high Pram also 
implies in a small Aeff [see Eq. (pT]) ]. 

The planet is supposed to be in the corotation radius, 
if the system is in a tidal equilibrium state (i.e., synchro- 
nized with the stellar rotation). However. iLevrard et al.l 
(|2009l ) have shown that in the majority of the observed 
transiting extra-solar planets, the star's rotation is not 
synchronous. Due to this reason, we compute r^^ for a 
range of equatorial radial distances. 

Considering the stellar winds from the simulations A, 
E, and F, we calculated at the equatorial plane by 
assuming that the giant planet has the same mass and 
radius as Jupiter, and that the planet's ma gnetic field is 
Bn = 100 G (same parameters adopted bv ILovelace et al.l 
l2008l ) . In order to perform the comparison of our results 
with the WD model, we evaluate Eq. ([12]) for both mod- 
els. We thus compute the velocity and magnetic field 
profiles of the stellar wind using the WD model. We as- 
sume the same parameters at the base of the wind as used 
in our simulations A, E, F, with the difference that Bq is 
the magnetic field intensity at the pole in our model, and 
in the WD model it is the radial field calculated at the 
surface. We do the same procedure for the more realistic 
winds we simulate, i.e., we use the results of the simula- 
tions for u, B, and p and evaluate using Eqs. ([TQ|) to 
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(fT2|) . These results are shown in Fig. [H where the sohd 
Unes represent the results obtained using our model and 
the dashed lines for the WD model. 

Figures [4^ and St) show, respectively, the effective 
cross-section of the planet and the ram pressure ex- 
erted on the planet for the region inside ~ 0.1 AU (i.e., 
< 10 i?*). It can be seen that when Pram is large, Aeff is 
small, and vice- versa. In addition, Aeff (Pram) calculated 
from our stellar wind model is larger (smaller) than if it 
were calculated using the WD model. This is because, 
for the WD model. Pram is verified to be essentially mag- 
netic, i.e., the first term of the right-hand side (RHS) of 
Eq. ([To]) is negligible compared to the second one, while 
for our stellar wind model, both terms are significant. 
In addition, for both models, the first term of the RHS 
of Eq. ([To]) are of the same order. The difference thus 
comes from the magnetic term of Eq. ([TO]) . that for the 
WD model is larger than for our model. This behavior 
is illustrated in Fig. [5l for case F. 

This results in larger time-scales as compared to the 
WD model (Fig. Hj^), which implies that the more re- 
alistic winds simulated in the present paper are not ex- 
pected to have significant influence on hot- Jupiter migra- 
tions. The time-scales we estimate (e.g., > 50 Myr for 
cases A and F, fo r r„ ^, 3 R^) ar e much larger than the 
ones estimated bv lLovelace et al.l (j2008f ) (-2-20 Myr). 
Other effects, such as an increase in the base density or 
magnetic field intensity, can reduce r^^. Case E presents 
the highest base density, while case F presents the low- 
est one among our simulations. From Fig. [4]3, we note 
that an increase in po reduces r^. We expect that for 
po > 2.4 X 10"-*^^ g cm~^ with the remaining parameters 
as adopted in cases A, E and F, the effects of the mi- 
gration will be more significant. A similar trend, i.e., a 
decrease in r^^, is also expected if one increases Bq. For 
example, cases A, E and F assume Bq = 1 kG. If such 
value were increased, it would result in a larger ram pres- 
sure exerted by the wind on the planet. This would cause 
a more efficient drag and ultimately it would lower r^^, 
and perhaps it would even become smaller than 20 Myr 
in the region where < 10 i?*. In this case, the more 
realistic winds could have a more significant infiuence on 
the giant planet migration, though it would still have 
much larger than the one compu te d wit h the WD model. 

According to iLovelace et al.l ()2008l ). the magnetic 
winds are expected to be important for planet migration 
for the case of a multipolar magnetic field rather than a 
dipole field. In the present paper, we extend this state- 
ment to the case of the multi-field component coronae, 
as the ones analyzed previously, where open and closed 
field lines coexist. Our model suggests that the stellar 
wind of these modified-dipole coronae are not expected 
to have significant infiuence on hot-Jupiters migration. 

If the wind azimuthal v elocity u^p is greater than the 
planet Keplerian velocity, ILovelace et al.l (|2008l ) showed 
that the stellar wind causes the planet to move outwards. 
In all the simulations we presented, this is the case only 
for case G, for > 1.26 For all the other cases, the 
wind would act as to push the planet inwards, and the 
closest the planet is from the star, the stronger /faster 
would be this dragging. 

An aspect that is not investigated in the present pa- 
per is when the surface dipole is tilted with respect to 
the stellar rotation axis. The wind ram pressure is an 
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Fig. 4. — Comparison of: a) The effective cross-section of the 
planet that intercepts the stellar wind, b) The wind ram pressure, 
c) The time-scale for the stellar wind drag change significantly the 
planet's orbit. We use the same parameters as in cases A (red 
lines), E (blue lines), F (green lines): solid lines are the results for 
our stellar wind model, and dashed lines for the WD model. 
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Fig. 5. — Relative importance of the first (dashed Hnes) and 
second (dot-dashed hnes) terms in the RHS of Eq. ()10|) calculated 
for case F using the results of our simulations (black) and using 
the WD model (red). Solid lines are the total ram pressures, i.e., 
the sum of first and second terms. 

interplay between magnetic (ex B'^) and kinetic terms 
((X p[U(^ — vk]'^). In the tilted case^ a planet orbiting the 
rotation equator would be subjected to high magnetic 
field intensities (since |B| is maximum at the magnetic 
pole and minimum at the magnetic equator). However, 
it is unclear the latitudinal dependence of u^: in the 
aligned case, u^^ is maximum at the equator; in the mis- 
aligned case and in a highly magnetized ambient, the 
stellar magnetic field is expected to channel the rotat- 
ing outflow. Hence, by a superposition of the magnetic 
and kinematic effects, it could be possible that be- 
comes smaller in a tilted configuration. A future work 
will explore this scenario. 

Another interesting configuration is when the orbital 
axis of the planet is inclined with respect to the spin 
axis of the star, as it is belie ved to occur, e.g., in the 
extra-solar planets HD80606b (Moutou et al.ll2QQ9D and 
X0-3b (Hebrard et al. 2008; Winn et al. 20091). In this 
configuration, the planet would interact with a stellar 
field strength that is larger at the position of the planet 
than at the magnetic equator. This could have an effect 
on that could reduce it. 

6. DISCUSSION AND CONCLUSION 

Recent works indicate the presence of complex mag- 
netic structure, as the one observed in the Sun, on T 
Tauri stars. This structure is certainly influencing the 
winds of these stars. WTTS offers a tool for studying the 
magnetic behavior of low-mass pre-main-sequence stars 
with the benefit that it is no longer strongly influenced 
by the presence of an accretion disk. As a consequence, 
the study of WTTS may eventually have applications on 
theories of CTTSs and also clues on the physical condi- 
tions for planetary formation and/or migration. 

In the present work, we showed 3D MHD numerical 
simulations of magnetized stellar winds in rotating T 
Tauri stars. With the parameters we selected, we scan 
the possible values of rotation period (0.5 d to 10 d), 
magnetic field intensities (200 G to 1 kG), and surface 
density (2.4 x 10~^^ to 2.4 x 10~^° g cm"^). We show 



that the plasma-/? parameter is a crucial parameter in 
the acceleration of the wind and in shaping the magnetic 
field lines. If the surface magnetic and thermal energy 
densities are equal at the pole, i.e. /3o = 1, we showed 
that: (i) the wind presents approximately spherical sym- 
metry; (ii) the magnetic field lines are not significantly 
distorted from the dipolar field we assume as initial con- 
dition of the simulation; (iii) the wind velocity does not 
differ considerably from the hydrodynamical case. On 
the other hand, if the magnetic energy density is much 
greater than the thermal one, i.e. Po <^ 1^ we observed 
departure of spherical symmetry. The wind shows two 
different components of the velocity, being faster at high 
latitudes and slower at low latitudes. In addition, the 
field lines are not anymore dipolar: the high velocity of 
the flow is able to stretch and open the field lines at 
high latitudes, while in low latitudes we observe elon- 
gated magnetic features. This indicates that streamer 
structures should be pre sent in WTTSs with 0o ^ 1 • 

In a different context, lud-Doula fc Owockil (|2002l ) de- 
fined a similar parameter that is used as an indicator of 
the effectiveness of magnetic fields in channeling outflow 
from hot supergiant stars. The "magnetic confinement 
parameter", defined by them as the ratio between the 
magnetic and kinetic wind energy densities, operates in 
a similar fashion as /3o operates in our work: it identi- 
fies whether the magnetic field will play a significant role 
in the acceleration of t he wind. Howe ver, both in the 
present work and in Vidotto et al. (2009), the thermal 
energy of the (coronal) wind is more significant then the 
kinetic energy, while in the line-driven winds studied by 
lud-Doula fc Owockil (|2002l ). the kinetic energy is more 
important than the thermal energy of the wind. 

We also show that when the self-consistent interaction 
of the field lines and the wind is taken into account, the 
picture of approximately radial streamers (near the star), 
as observed in the Sun, is not anymore verified for the 
WTTSs with rotation periods up to 10 d: both open and 
closed lines are twisted. 

With typical parameters of WTTSs, we estimate the 
time-scale for planet migration due to the torque ex- 
erted by the azimuthal ram pressure of the stellar winds, 
and we showed that for the multi-component corona of a 
WTTS simulated here, this migration m echanism is not 
important, in contrast to the results from lLovelace et"al] 
([20081) , who argued that magnetic winds computed using 
the WD model are expected to be important for planet 
migration. Further simulations with a greater parame- 
ter range (higher magnetic field intensity or higher den- 
sity at the coronal base) or different system configura- 
tions (a tilted stellar magnetosphere or a tilted plan- 
etary orbit with respect to the stellar spin axis) may 
reduce the time-scale to significantly affect the orbital 
radius of the planet. For instance, if we consider a mag- 
netic field intensity of 2 - 3 kG, still consistent with 
values derived fro m Zeeman-broadening measurements 
() Johns-Krullll2007l ) . a reduction in is expected, as the 
ram pressure of the wind increases. 
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